function r=f(t, y)
    TC = 0.1;
    TV = 0.2;
    r = TC * t - TV * t;
